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I. INTRODUCTION 


A nuclear pumped laser, operating at a wavelength of 

3 

1.79 pm. on the 3d{l/2)-4p(3/2) transition in argon with He 
as the majority gas has been produced by dalufka et. a1, at 
NASA-Langley Research center. A model, proposed by Harries 
and V/ilson, to explain the operation of this laser uses col- 
lisional radiative recombination of Ar as the population 
mechanism of the 3d{l/2) state. This model also proposes 
dissociative recombination of Ar^. as a population mechanism 
•for the 4p(3/2) state. 

■ Work, done by Biondi et. al . at the University of Pitt 
sburgh. confirms that' dissociative recombination populates the 
4p level but indicates that dissociative recombination domin- 
ates oyer collisipnal radiative recombination in pure argon 

■3 

at low pressures. As the laser developed at NASA-langley 
■operates in the high-pressure regime, calculations have been 
perTormed by DePaola et. al . (see Appendix 1) extrapolating 
v/£Tl-known low-pressure results to the hi gh-pressure regime. 
From the results of these calculations there is some question 

-in predicting- laser operation 'at T.79 pm. in' a ^He-Ar plas-ma^. 


The purpose of this study is /the investigation of the 



energy pathvrays in He-Ar gas mixtures at lov; pressure is acc- 
omplished by observing the effects of varying partial pressures 
on the emissions of levels lying above the 4p level in argon 
during a pulsed afterglow. An attempt v/as made to investigate 
the population mechanisms of the 3d level in pure argon by ob- 
serving emission from the 3d ( 1/2 ) -4p (3/2} transition in a high- 
pressure plasma ejccited by a high-energy electron beam. 
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II. THEORY 

• 3 

A nuclear pumped, He-Ar laser, excited by the 

He(_n,p}"H reaction, was created by Jalufka et. al . in 1976. 

The laser operates at a wavelength of 1.79 pm. on the 3d(l/2)~ 

4p(3/2) transition in atomic argon. A model for the laser, 

proposed by Harries and Wilson, assumes that the population 

mechanism of th.e 3d(_l/2) state is collisional radiative recom- 

bination of Ar and that dissociative recombination of Ar^ 

populates the 4pC3/2) state. ^ 

A flow diagram, illustrating the energy pathways 

proposed, ts shown in figure 1. The greater part of the fis- 
sion energy ' produces helium ions. He , and excited helium, 
atoms In a, n)etasta.hle state, He(2^S). Some of the He(2^S) is 
converted to a molecular helium metastable, He2(2 2). Argon 
ions are created from both atomic and molecular helium meta- 
stables by Penning ioni 2 ation. The helium ions are . convert- 
ed to molecular helium ions, He^, which also serve to ionize 
the argon by cha.rge transfer collisions. The- ionic argon 
may recombine by collisional radiative recombination or be 
converted to Ar^. The products of collisional radiative re- 
combination cascade through excited states and may produce a 
population inversion between the 3d(l/2) and 4p(3/2) states. 
Molecular argon ions undergo dissociative recombination which 
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Figure 1. Flow Diagram of Processes Important' to the He-Ar 
Laser. 

Dissociative recombination is assumed to populate 


the lower laser level. 















is presumed to populate the 4p(3/2) state thus becoming a sr 
ource of population for the. lov/er laser level. 

The assumption that dissociative recombination popu- 
lates the lower laser level is supported by the work of Biondi 
ans Shiu v;hich shows that the ■ip states are the most strongly 
populated. Experiments performed by Frommhold and Biondi, 

•in '1 969, indicate that dissociative recombination dominates 

over collisional radiative recombination in pure argon ttilow 
4 

pressure. 

Calculations performed by DePaola et. al. (see Appen- 

O 

dix 1) provide a characterization of high-pressure He-Ar 

plasmas based on extrapolation of the knowledge of similar 

low-pressure plasmas? Results of these calculations for a 
3 

He-Ar plasma at a total pressure of 2 atmospheres with a 1% 

1 8 

concentration of a'rgon and a thermal neutron flux of 10 neu- 
- 2 - 1 

trons-cm -sec are quoted . here 

,Gas Temperature T^=300 K 

Electron Temperature T^=536 K 

14 -3 

Electron Density N^=2.'08X10 cm 

kr^ Density : {Ar‘^}=7. 57X10^^ cm'^ 

Ar^ Density {Ar^l^l . 08X10^^ cm 

Collisional radiative rate coefficients for the minority gases 
in various helium-dominated gas. mixtures have been calcuated 
by Whitten et. al . ^ For a He-Ar mixture at 2 atmospheres 
total pressure the rate- coef f ici ent is: 

55X10"^ cm‘^^-sec“^ 

• * 

and the rate coefficient for dissociative recombination of 
is : 

a-p=9.1X10~^(T^/T^)“- cm.^-sec"? 


^ POOR QuALrry 
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The assumption may be made that all of the products of coll- 
isional radiative recombination pass through the 3d level and 
are distributed according to the densities- pf the i nd i v i dua 1 ist- 
ates'within’the level. The same assumption is made with res- 
pect to the products of dissociative recombination cascading 
through the 4p level. The rates of -poulation of the two 
states are thus: 

d{ArC3d{l/2))}/dt=2/20Xa^j^j^{Ar'^}Hg 

=4.02X10^^ cm“^-sec"^ 
d{Ar{4p{3/2))}/dt=4/12Xap{Ar2}i'-!g 

=4.78X10^^ cm"^-sec"^ 

The rate of population of the lower laser level is thus seen 
to be three orders, pf njagnitude greater than that of the upper 
laser level and laser action is not predicted in this model. 

The character i.sti cs of an existing He-Ar laser may 
be used to predict the value of the collisional radiative re- 
combination rate coefficient necessary to produce lasing. . A 
laser described. -by DeYoung et. al. operates at 2 atmospheres 
total’ -pressure ..v;ith a 1 % argon concentration. • The peak pov/er 

output of the laser is 3.4 Watts for a thermal neutron flux 

17 -2-1 ■ 

of 10 . neutrons-cm -sec . The efficiency of the. laser 

is reported to be 0.09%. The power depos.ited in a laser of 

this type is: 

where a=5X10 cm is the fission cross-section of He, ^ 
is the thermal' neutron flux, Ep=0.79 Mev is the energy re- 
leased per fission event and V is the active volume. The 



power output is given by: 


= ea<|){^He}EpV 

w,here e is the efficiency of the laser. Thus the active vol- 
ume is given by: 

V=r^/eaM^He}Ep 
= 1.162 ctn^ ' 

1 9 

Th.e peak power output of 3.4 Watts corresponds to 3.06X10 
- 1 3 

photons-sec at 1.79 pn)- As the . He-Ar laser is a steady- 
state laser the ^source term for the population of the 3d(l/2) 
state must be equal to the loss term. Using the model of 
Harries and Wilson this requirement becomes: 

cccRRNg'{Ar’^}V=3. 06X10^^ 

or : 

3 . 0 6 X 1 0^ N g { Ar '*■ } V 

Th.e computer program used by DePaola et. al . was modified 
slightly .Csee Appendix ’2) to calculate the electron and ion 

densities of a Ke~Ar plasm at 2 atmospheres with a thermal 

■ >7 ? - - L ■ ' 

neutro.n flux., of 1.0 neutrons-cm -sec . The results are: 

Electron Density N =5. 53X1 0^^^ cm"^ 

Ar^ Density {Ar+}=1 .’SIXIO^^ cm“^ 

And the required coll isional. radiative rate coefficient is: 

63X10'® 

a full order of magnitude greater than that calculated by 
Whitten et. al.® ' 



9 


III. EXPERIMENTAL PROCEDURES AND EQUIPMENT 

The nn'xtures of interest in this study v/ere excit- 

ed by a pulsed electrical discharge. Pulse duration v/as 4.5 
±.05 ysec with a peak current of 40+5 amps (see figure 7). 
Electrode spacing v/as 7.5 cm in a glass discharge tube with 
quartz windov;s. The discharge tube was 2.5 cm in diameter and 
1 m long 


figure 2 shows a block diagram of the vacuum and gas- 
handling system. The pumping system consists of a 4-inch oil 
diffusion pump with a liquid nitrogen cold trap and auxiliary 

roughing pump. The diffusion pump was capable of bringing the 

- 7 

entire gas-handling system .to a pressure of 10 Torr and the 

_ 5 

discharge tube vfas brought to a pressure of 10 Torr or 
better prior to each gas fill. | 

The gas inlet system consisted of reservoirs of re- 
search grade helium and argon with control 1 ed-fl ow valves conr 
nocted to the remainder of the system through a length of 1/4- 
inch copper tubing. Gas pressure was monitored by an absolute 











pressure gauge calibrated from 0.1 to 20 Torr in 0.1 Torr gra- 
dations and a differential gauge cal i brated from 0 to 800 
Torr in 10 Torr gradations attached to the gas inlet system. 
Under vacuum conditions pressure could be monitored by the 
mass spectrometer or by an ionization guage (Varian model 563) 
calibrated from 10'*^ to 10"^ and attached to the gas chamber. 

Preparation of Gas Mixtures 
In preparing th.e gas mixtures to be studied the en- 
tire system v/as first flushed v/ith the minority gas to re- 
duce impurity levels. The minority gas was used in flushing 
as insurance against possible incomplete pumping due to the 
lov-^ conductance of the 1/4-inch tubing used in the gas inlet 
system. Following this the discharge tube was. overpressured 
v/fth the- desired gas. mixture and subsequently pumped to 10. 
Torr. Before each run the. discharge tube was isolated from' 
the vacuum • and' gas-handl i.ng system. 

Data Aquisition System 

The data aquisition system is shown schematically in 
figures... Data v/as accumulated by a photomultiplier tube (RCA 
■model 7625) attached to a 0.5-meter Ebert scanning spectrome- 
Iter . (Oarrel -Ash model 82-000). The spectrometer was calibrat- 
ed to within 5 angstroms with a mercury arc lamp and location 
of individual lines v/as verified by maximizing the signal 




Figure ,3. Block Diagram of ..Data Aquisition System. 

The signal from the PMT v;as routed through a pho 
-ton Di scrimi nator/Ampl 1 fi er before goi ng to the 
photon counter. 







■from the photomultiplier as viewed on an oscilloscope. The 
photomultiplier signal was then routed to a photon discrimin- 
ator/amplifier- (Mech“Troni c s model 511} to a photon counting 
system. 

The photon counting system consisted of a photon coun 
ter (Ortec model 9315) with a sampling/control unit (Ortec 
model 9320). The photon counting system was operated in a 30 
step boxcar mode allowing counting to take place during a pre 
selected time period at 30 different temporal positions in 
the active discharge and afterglow (see figure 4). A pre- 
trigger was supplied to the sampling/control unit by the dis- 
charge pov/er supply trigger 1 msec before the initiation of 
the current pulse. A variable delay line in .the sampling/- 
con'trol unit served to initiate counting approximately 20 
psec in advance of the current pulse for each cycle. During 
step.l the photon counter is enabled by a 40 ysec gate sig- 
nal supplied by the sampling/control unit. This is repeat- 
ed foreach trigger pulse received by the sampl ing/control 
unit for the duration of an- .integration ’interval of 5000:'c'urr 
ent .pulses. At the end of the integration interval the data 
accumulated by the photon counter is transferred to .a buffer 
where it is available for recording. The delay to the gate 
signal is then incremented by twice the selected gate width 
and step 2 begins. Operation during step 2 is identical to - 
that of step 1 except that counting begins at a later time 
in the discharge/afterglow. A gate width of 40 psec gave 
samples as far as 2.4 msec into the afterglow. An integ- 
ration interval of 200 sec yielded over 5000 current pulses 




Calibration of the Optical System 

The spectral response of the spectrometer and photo- 
multiplier tube v/as calibrated v/ith a tungsten filament lamp 
(ITT model T-lO) placed behind the discharge tube, the temp- 
erature of the filament v/as 1500 K. The spectrum of the lamp, 
corrected for the self-absorption of the filament, is showm 
in figure Sand the spectral response of the optical system 
is shown in figure 6. The response curve shows the signal 
■from the photomultiplier tube, as a function of wavelenth, 
relative to the response at 3600 angstroms. 

Electron Beam System 

The system used for high-pressure measurements is 
shown schematically in figure 7. It consists of an electron 
beam machine supplied' by Maxwell Research Co. attached to a 
stainless steel plasma cell 2S..5 Inches long and 12 inches- in 
diameter.. Twa CaF 2 - wi ndows , 2. inches in diameter, were att- 
ached to the cell at its midpoint. The electron beam deliver- 
ed' approximately 7.5 amps-cm ^ at 60 kV in 900 nsec pulses. 

The gas-handling system .and preparation procedures were sim- 
ilar to those described above v/ith the exception of the to- 
tal pressures. A pressure of 300 Torr v/as used for the elec- 
tron beam measurements. Data v/as accumulated by a liquid nrt- 
rogen-cooled InAs infrared detector (Barnes Engineering model 
0211) attached to the same spectrometer described above. The 
spectrometer v/as equipped with an infrared grating (Oarrel- 
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Figure 5 Spectral Response of Calibration Lamp. 

Blackbody curve of a tungsten filament at 150C 
corrected for self-absorption. Intensity is ' 
arbitrary units. 



INTENSITY 


SPECTRAL RESPONSE OF CALIBRATION LAMP 
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Figure 6 Spectral Response of the Optical System. 

Response relative to the response at 3600 angst- 


roms is shown. 





SPE-CTRAL RESPONSE OF OPTICAL SYSTEM 









Ash S1347D8A9). The signal from the detector was ampilified 
by a 60-db gain low-noise preamplifier (Ortec model 7506) and 
viewed directly on an oscilloscope (Tektronix model 466) in 
storage mode. 



IV. EXPERIMENTAL RESULTS 


The results of this stud.y are presented in three sec- 
tions. First, the intensities of observed emission lines are 
(presented. A quantitative discussion of the results is pre- 
sented in the s.econd section. This is follov^ed by a summary 
of the conclusions drawn from this work. 

Intensities of Some emission Lines in Argon 
Emissions of lines emanating from the 7s , 6p , 6s , 5d ,5p , 
4d and 4p levels in argon v/ere monitored during the active dis 
charge and afterglow in three different gas mixtures. The gas 
mixtures studied v/ere pure argon,, hell-uni with 1% argon and 
•a 50-50 mixture of helium and argon. Plots of photons acc- 
umulated vs. time are shown for several, representative lines 
in figures 8-10 The relative intensities of each of the 
lines are tabulated in tables 1-3. Relative intensities are 
reported for the active discharge and 5 times during the after 
glow. Transitions chosen for analysis are the 4p'(l/2)- 
4s-’ (1/2) transition at 7503 angstroms and the 5p'(l/2)- 
4s’ (1/2) transition at 4259 angstoms. The 7503 angstrom line 
was chosen because the 4p' level lies well below the ground 

a. « 3 

state of Ar 2 and is one of the levels' of interest in the He-^ 
Ar laser. The 4259 angstrom line was chosen because the 5p' 
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Figure 8. Plot of Intensity Versus Time at 7503 Angstroms. 
Uncertainty is ±0.14X(Intensity)y“ 

Symbols are: +-pure argon 

x-1% argon 
*-50% argon 

The current pulse occurs at' 20 psec. 
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Figure 9. Plot of Intensity Versus Time for 58S8 Angstroms- 
Uncertainty is ±0.14X(Intensity)^(^ 

Symbols are: +-pure argon 

X-1% argon 
*-50?S argon 

The current pulse occurs at 20 psec. 
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Figure 10, Plot of Intensity Versus Time for 4259 Angsirrons^ 

Uncertainty is t0.14X(Intensity)^{^ 

Symbols are: +-pure argon 

X- 1% argon • 
argon 

The current pulse occurs at 20 psec. 
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Time after trig cjor (psec ) 


Ua vel ength 
(angstroms) 

0 

80 

1 60 

240 

320 

400 

7503 

233.91 

83147 

71344 

58626 

2739 

6268 

7030 

12352 

28019 

1 3232 

3709 ■ 

4S2 

47 

6937 

3712 

1 5929 

5003 

2386 

761 

T88 

5888 

4431 

5710 

2386 

662 

1 23 

18 

:5650 

449 

3869 

761 

395 

97 

17 

'^i259 

1382 

2348 

188' 

365 

72 

- 1 1 

3649 

698 

93 


11 

3 

1 


Table 1. Relative Intensities in Pure Argon. 

Trigger occirrs 20' ysec before current pulse. 





Time 

after trigger 

(nsec) 


Wavelength 
(angstroms ) 

0 

80 

160 

240 

320 

400 

7503 

12885 

174545 

1 5743 

584 

31 

37 

7030 

42236 

35455 

2143 

129 

25 

14 

6937 

6612 

22857 

1688 

128 

29 

21 

5888 

10590 

5570 

333 

49 

30 

22 

5650 

91 56 

3308 

18 

. 19 

8 

3 

4250 

1158 

4887 

321 

17 

4 

2 

3649 

1630 

100 

10 

3 

2 

1 


iable 2. Relative Intensities in 99% Helium With 1% Argon. 


Trigger occurs 20 ysec before current pulse. 
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Tinie after tricjger (iisec) 


Wavel ength 
(angstroms ) 

0 

80 

160 

240 

320 

400 

7503 

17588 

367639 

303431 

158741 

27582 

3535 

7030 

6719 

1 20476 

41412 

4429 

264 

33 

6937 

10184 

63237 

4362 

415 

60 

55 

5388 

6285 

1057 

241 

86 

35 

18 

5-650 

9637 

1397 

397 

64 

16 

7 

4259 

982 

10672 

3713 

574 

61 

10 

3649 

220 

142 

41 

12 

5 

1 

Table 3. Relative 

Intensities in 50!^ 

Hel i um 

w i t h 

50% Argon 


Trigger occurs 20 ysec before current pulse. 
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level lies near the 3d level. 

The evolution of the chosen' transitions may be examined 
through the ratios of their emissions. The ratio of emission 
at 4259 angstroms to that at 7503 angstroms in pure argon falls 
off- from 0.059 during the discharge to 0.028 50 psec after the 
discharge andO.005 130 ]jsec after the discharge. A similar 
situation is encountered in the 1^ and 50/5 mixtures. This is 
to be expected. 

As the partial pressures of the gas mixtures are al- 
tered a marked effect is observed. The ratio of emission at 
4259 angstroms to that at 7503 angstroms rises from 0.059 in 
pure argon to 0.090 in helim vn'th T% argon then falls off to 
Q.05G in the 50" mixture^ during the'.active discharge. Early 
in the afterglow the ratio remains nearly constant as the par- 
tial pressures change but by 130 ysec into the afterglow the 
behavior is similar to that of the active discharge. 

The high-pressure experiment yielded indeterminant re- 
sults. Infrared radiation from an electron beam-generated 
plasma of pure argon at‘300 Torr was scanned over the range of 
T . 5 ym to 1 . 8;4 ym. No signal was observed above the radio- 
frequency noise generated by the electron beam machine in- - 
dicating that the sensitivity of the I.R detection apparatus 
was insufficient. 


Discussion 

The plasma is assumed to be in a quasi-steady state, 
thus the intensity of emission from a given state, as expressed 
in number of photons is equal to the rate of population of 
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that state. In the model of Harries and ’Jilson dissociative 
recombination is assumed to be a population mechanism for the 

p 

4p level but not for the 5p level. CoHisional radiative re- 
combination is assumed to be a population mechanism for both 
1 evel s . 


The collisional radiative recombination rate coeff- 
icient is typically expressed as the sum of a neutral part and 
an electronic part. In the 1 ov;-press ure regime it may.be ass- 
umed that the total rate coeffiecient is approximately equal 
to the electronic parf^ The collisional radiative rate coeff- 
icient is expressed as': 


where A, 


^CRR^^^CRR^’e^e 

is a constant and x is unknov/n. The dissociative 


recombination rate coefficient is 


— A T ” • 

Thus the ratio, R, of emission at 4259 angstroms to that at 

I 

7503 angstroms' may be expressed as: 

The quasi-steady state assumption imposes the requ 1 rement- that: 
d{Ar2>/cit=k{He}{Ar}{Ar'^}-AQT^^'®^N^{Ar2} = 0 
where k{He}{Ar}{Ar } is the rate of formation of 

{Ar‘^} = k{He}{Ar}{Ar2} 

A T 

^D'e '^e 


Thus : 


Aj,ppTe“^Ne^{Ar‘^} + k{He}{Ar}{Ar'^} 
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R= T 

Uk{HO{Ar}/A(.p,^T/>‘N^2 

’’'“ ^CRR^e ’‘**6 
k{He}(Ar} 

In this expression {He} and {Ar} are constant uith 

respect to time.- Thus the ratio of the ratios at two differ- 
ent times during the afterglov; may be written: 

R-j/R2=CNg-j/Ng2^^^T’e}/^e2^ 

and : 


-x = ln{(R^/il2)(Ng2/S'*el5 

The time evolution of ’electron density and electron temperat- 
ure in a plasma under conditions similar to those in this study 
v;ere measured by Monchicourt et. al.^ Their results are shown 
in table 4 along with the other data pertinent to this dis- 
cussion. The calculated values of -x are found to be 4.05 for 
pure argon, 0.255 for argon and 6.43 for 50i4 argon. 

Concl usi ons 

Results of the high-pressure observations are incon- 
clusive, -indicating a need for further work in this area. The 
electron temperature dependences calculated from the results 
of the 1 ovz-pressure observations, on the basis of the model of 
Harries and Wilson, are incompatible with generally accepted 
values. The model of- Harries and Wilson and the results of 
this study are seen to be inconsistent v/ith each other. 


ORIGINAL PAGE IS 
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Time after discharge 

50 ysec 

130 psec 

R in pure argon 

0.028 

0.014 

■ R in 9 9 % He v; i t h 1 % A r 

0.028 

0.020 

R in 50f. He with 50 % Ar 

0.030 

0.012 

g 

Electron Temperature 

780 K 

710 K 

8 

’El>e.ctron density 

9.0X10^^ cm“^ 

X 

1 .^ 

o 


Table 4. Parameters Used in Calculation of Tempera ture De- 
pendence of Col 1 1 si onal Radiative Recombination 
Rate Coefficient. 



APPENDIX 1 


ESTIMATIONS 


D. DePaola, 


OF ELECTRON DENSITIES AND TEMPERATURES 
IN ^HE-DOMiriATED PLASMAS 
by 

S. D. Marcum, H. K, Wrench, B. L. l/hitten 


and W. E. Wells 
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•ESTKiAiIC.'iS OF ELECTRO.'! OETSITIES /TlO TEM?=.RATL'?ES 
3N ^He-00:ii:;AT£D PLAS.'IAS* 
by 

8. D. De Paola, S. D. Marcuirij K. V.’re.>:ch, E. L. llhittS.T 
- and IJ. E. Kells 


IflTROD’dCTICM 

It is veil recosnl2ed that in order to make efficie.’it use of the large 
-energies derived fro.n a fission reactor ‘..’ith a iiuclear PuT-aacl Laser (URL), a 
medium that v/ill lase at high pressure must be en'plcyed. This is only one 
of the o-any features of future PiPL’s that will make them oarksdly different 
from coiTJEon lasers of today. ✓ With this in mind ue have Bade an attei^pb to 

3 

characterize a high-pressure-, He-dcminatad RPL plasma, based on what u kno.>n 
or similar low pressure plasmas. 

i’lany processes important to i.'PL’s are strongly depe.adsnt upon electron 
temperature and/or electron densit3*. Also, due to the nature of nuclear- 
pumping, performing diag.i-ostics on fJPL plasmas is a very difficult, expensive 
and hazardous task. These facts provide the impetus for perfornir.j the calcu- 
lations described here. These calculations are straightforv/ard extrapolations 
of lew pressure results into the high pressur-e regime. Several of the uncer- 
tainties inherent in this approach ere a1.so discussed. 

THE PUP.Z'^He PLAp^. -■ 

Taking helium as the'sole species in the -plasma, the following reactions • 
\.‘ere considered: ■ 


1) S = c; n 


C^He)^ 


Eouaticn 1 describes the production or He whe.-e o is "tr.e lissicn cro3j.-sect'vn 


for is the thermal neutron flux (asscTed to be 10^^ n-cia '^-sec ‘), '’s 

the densiti' of .is, e- is the cnergj' released per fission of Ke and If is the 


-2 
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, energy investoc! p2>’ clectroti-ion p<iir proiiac-.-d (=-'.T.5eV). 

2) = S/.56. 

This equation accounts for production of s) r-tastablss which 
proceeds at roughly twice the rate of Ke production. 

3) He"^ + 2He — > iie^ + he. 

3~faody conversion to He 2 > This process varies as 72^^,^ (average of 
several fnaasurements). 

4) Ke(2^S) + He(2^S) — ^ Ka"^ - Ke + e. 
Kstastable-fr.etastable ionization which has the rate coefficient 

B"1.8xl0^ ctn^-sec 

5) Ke(2^S) + 2He — •> He2(2^Z) + He. 

' 2 
3-boriy conversion of trstastahiss i;h1ch scales as 

6) -J- e +'^ CRR. 

This represents col lisional -radiative recombination which includes a 
very large number of passible processes. 

7) Hs(2^S) + e — ? He + e. 


The last process considered is superelastic relaxsticn. The rate for this 
reaction is temperature dependant (T^ ' } and \;as derived from £ detailed 
balancing calculation based on the forward process. (Reactions involving the 
helium singlet iratastable were not considared because being only .8eV away 
from the Ha(2 S) state 'it is very, quickly converted to the triplet Ketastabla). 


Using the above processes the following rate equations were developed: 

• S>l/2 p.[He;2%)] - 78 [Ks^j 

■Hare’ the" first' source- term is due to ionization 'by high' energy electrons 
produced -by' tha, fission products of-^Ha; the -second is due to retas-table-ratastabls 

-f. 

ionization. The loss term is due to-3-body conversion. For He 2 we hava 


cjrtiet] 


~ = 73 P„^ [He 3 - a [Ke,] n , where 




ere electron and cas temperatures and tn- electron end rieutre’ densities, ' 
respectively. The source tonn is due to 3-body conversion v.-.hereas the loss is 
due to collisional-rediative recosibinaticn. h'ota the electron ten.perature end 
density dependencies. For heliira n’etastables the rate equation is 

dCHeil^] = ^ ^ _ g - Y [:-!e(2’s)] - .6 [}le(2^S}3 

uhara v = 7 X 10 -j-V2 _ source terns here are due to excitation by 

e 

3 

high energy electrons produced by the fission products of He and collisional- 
radiative recombination of Ke^- Losses are due to ir-stasteble-rratastable ionization, 
superelastic relaxation and 3-body conversion of retastables, respectively. For 


electron density v;e have 


i!!§ 

dt 




=. S + 1/2 g [Ha(2^S)j- - o [Hep n 


•Here, of course, vie have the sarre source terits as with Ha ; the loss is assu.rad 

to be due to collisional-radiative recombination. Also 'considered was the 

following energy balance equation; 

d(3/2 nkTl = S -h 1/2 0 [He(2^Sl- A='e_^ ^ Y 

dt 

^(3/2 t) (T^ - T„) 


The source terms are due to ionization by high energy electrons produced by the 
3 

He fission products, iretastable-metaotable ionization and superelastic 

relaxation- As A s‘ , and A Ep ere defi.ned as the energies transferred 

to the traxwellian part of the electron energy distribution by e non-maxwell ian 

electron at 10, 15 and 20 eV respectively. The loss term is due to electron' . 

neutral collisions with w being the electron-neutral collision frequency. 

eo 

We have an electron energy distribution like that of Figure 1. High 
energy electrons produced in REGIOil III by the ^Ke fission products cascade 
dow.n in energy toward a primarily pax'./ellian region of the distribution (RESIOm I)- 
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Bsyond the first excitation potential of He, inelastic collisions donJinate. At 
and below the first excitation potential of helium production of lie"'’ and He{2^S) 
occurs, '..'ells et. al. have shown that the energy transferred to the irax‘..’sl1 ian 
part of the electron energy distribution by a ncn-rax>..allian electron is a function 
of electron density and pressure [Figure 2). All values of i s used in this 
calculatio.n v/ere extrapolated from the work of 'dells. 

Using the above rate equations in steady-state form '..e calculated electron 
temperature and density as a function of pressure from IGO to 2330 torr. With 
these temperatures and densititas v»e \;ara also able to calculate the densities 
o'f Ha , Hfi 2 and the density of the Ha[2 S) metastable, all as functions of 
■pressure. However, pure heljum is not expected to ma’se a very good lasing rediuTti, 
so we next considered a mixture of helium and arcort- Results for the pure 
helium plasma ^Jill be shown later in comparison to those for the helium-argon 
mixture. 


THE Ka-Ar (li'Ar) PlASHA 

Upon adding 1" argon to the plasma, all reactions considered earlier with 
helium alone are still operative. Ho’.jever, several other reactions with argon 
must now be taken into account. In Figure 3, equations 1-7 are for a p'ure 
helium plasma. Equations 8-15 include possible He-.Ar reactions as 'well as Ar-Ar 
reactions. Eq. S - charge exchange with He^ ; 9- charge exchange with He^; 10 - 
3-body charge- sxchencs; T1 - Panning ionization of argon; 12‘- dissociative 
recombination; :'.T3 - 3 body- conversion;' H - col lisional-radi stive racembi nation; 
and 15 - associative Ionization to In the case of collisional-radiative 

recombination we assuced this reaction to be idantical to collisional-radiative 
recombi natl on of He. 

The rate equations listed earlier must now be modified to include these 
reactions. Figure 4 lists the updated set of rate equ-ations; additions to the 
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rr.te cq:.ations for the pure he'Sii.n plesira ere ii'Vj’erTir.a-J, For ozenolej Ke^ 
can no.: be lost through charge e^ichange uith Ar. he^ can also be lost thro::gh 
charge exchange as ‘..’ell as 3-hody ccnyersion, etc. However, considering steady 
state conditions , those ecjuaticns can be solved for electron terperature and 
density as functions of pressure, as well as the densities of all other species 
in the plasma. 

results 

Figure 5 is a plot that compares electron tenperat-jre with and ',/ithaut argon. 
Both are plotted against pressure from 100 to 23C0 torr. t!ote that the electron 
temperature is over ICCOK at low pressures, but tends toward roam tesrperaturs at 
'Mcher pressures. This is important because of the cfape.ndencs of the various 
coefficients on electron temperature. Also note that the addition of 1Z f\r to 
the plasma does not significantly affect the electron te:rperature. 

Figure 6 conmares electron densities with end ’.rithout argon. Here the 
variation in electron denisty uith pressure is not so great as was the temperature 
change, but again note the dependence of the rate coefficients on elactron density. 
Figure 7 is a plot of He"^ density versus pressure with Od end IS Ar. Upo.n adding 
argon to the plasma the He' density is lower at all pressures. This is rsesonabls 
considering that charge exchange between He and Ar can now occur. Similarly, ’ 
Figure 3 shews r.S;^ density with and without argon plotted against pressure. The 

• J. 

-large decrease in K=^ concentratlcn at all pressures is due to 2- and 3-bcoy 
charge exchange with argon. Figure 9 plots the helium metastable (f;e(2 S)) density 
versus pressure also with and without argon- The decrease in the helium meta- 
stable density with li; Ar at all pressures is due to Penning ionization of argon. 
filso, densities of Ar~, Ar^, and Ar '.;sre calculated as functions of pressure end 
are s;;cv;n in Figure 10- The argon excited state density represents an acccKulation 
of all ergon excited states that exist between the ionization limits of the atomic 
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ind inolscular ions into a single state. This v.'as dsne so that associative 
ionization of argon could be considered. 

raUCLUSIOilS 

The above results are reasonable based on extrapolation of tha knowledge 

3 

of lo-.;-pre5Sure Ke-dominated plasira s. However: 

1. What tha callisicnal-radiati vs reccr-binatiDn rata coafficiant at 
high pressures? V.'e have eraployed a low prssssura coefficient in our calculations 
that may or may not ba valid at higher pressures. 

2. We have assumed that reccmbinaticn of Ar^ is in the collisional 
radiative fomr, which seems reasonable, butvie’ve assumed collisional radiative 
recombination of argon to ba‘ identical to that of helium. This later assumption 
may well be a dubious one. 

.3. What other processes are important at higher pressures that we could 
not consider? It is possible that Ke^, Ari end/or He-r achieve large densities 

%j -a 

at higher pressures such that they become important species in the plasmas. These 
could not be considered in our analysis. 

These are exaiaplss of the questions that this type of calculation raises. 
While it is true that there are nsore sophisticated retnods of calculating electron 
temperatures end densitiss» given the uncertainties in cross sections and other 
data \;e feel that the approach outlined here is s'ufficient for the presant- 


*This research v.'£S s?cinEa.''ed in part by tha National Aeronautics and Space 
Administration. 
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REr.crioris cgnsidep.ed 


’• • (Hsnc„„e) 

2- = 5/. 55 (HaC2^S) source} 

3. Ke^ ZHe •> Ke| -i- He 

4. He(2^S) -5- Ha{2^S) Ha^ -J- He -!• e 

5- He{2^5) + ZHe > V.3^{^Z) ^ Ke 

6. K^-i- e i- X = 1 ^ CRR 

7. HeCZ^S) e ^ He + e 


-. p = 1.8 X 10^^ c(n^-5ec“^ 


« . 6 ? 


He 


7 X 10~^^ 




ETHER REACTIONS WITH 

AR5CW n-i) 

S. He^ A> > Ar- Ka 

V= 10-17 cm2 

9. Ke^ •*■ Ar ■ — > 

Ar”^ ^ 2H& 

pj= 2.2 X 10~1° cKi^-sec"! 

10. Ke^ -J- Ar -5- He 

. — K Ar^ -!- 3Ke 

= 2.4 X 10“^^ cni®-sac“^ , 

11. Ke(Z^S) + Ar - 

— j Ar + He + e 

&3 = 7 X ia"l® a/ 

12. Ar^ + B ^ 

d- 

Ar + Ar 

^=9.1 X iF {3G0/fg}*®l cn^-ssc-l 

13. Ar^ -e Ar + Hs 

V ArJ + Ke 

- 37S?HeV 

14. Ar^ -J- e -f- X : 

CRR 



2 X JO ca^ -sec~^ 


Figure 3: Reactions Considered in this Study. 
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RATE EQUATIO.NS 

= S 1/2 p Die(2^S)3 - 7S [He"] " 


dt 


= 78 [He"^] - a [Hsjl n^ - [ Hs|] [Ar] - [He|] [Ar] [He] 


. 78 [He-"] - - ^ 

2 

fe( . 2 S)] ^ 5^ ^ ^ p [He{2-"5)]" - y [He(2^S)] n^ - ,$ [Ke{2^S)] 


dt 


'» 3 ~\- i 2 


- ^3 [Ar] lHa[Z^S)2 


dt 


* * ^ 
= S -M/2 p tHeCe^S)]^ - . a [HeJ] %- [Ar] [H5(2 5)] - [Ar^] 


- c Ea/] n^ i 2 X 10“^ EAr’^] [Ar] 


o<v> [Kg'"] [Ar] -h EHe^] [Ar] Er:e|] [Ar] [He] -i- P3 [He{2^5)] .[Ar] 
- 375 Pj,^ P^^ [Ar^] - G [Ar ] • 

4[A)~o] q * ■ A • 4 

^ = 2 X 10 ^ [Ar ] [Ar] •^ 376 P^„ [Ar ] - [Arp 


dt 


3 - 2 X 10“® [Ar*] [Ar] - [Ar*]/a 


Khere: 


S = cr ^ Sc 

S^-S/.ES; ■ 

T = 7 X 10~-^1^ j]/^- ■ 


a<v> = G X 10 


i—20' f~r- i-r \ 5 

e{\^ «e 
13 


G = 4:5 X 10“'^^ CT„/T„) -*X 5 X 10-"^(T /Tj”'-’n 


-27,, j-l±l 
e o-' o 
.-11 


= 2.2 X 10 


,-10 


Pg “ 4.2 X 10 
Cg = £.1 X 10“^ (VT^)-®"' 


y. = 2-4 X 10' 


,-24' 


Figure 4: Rate Equations Employed in this Study. 
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Figure 5 : Plo!: of electron Te^iperatuce versus Pressure. 


ELECTRON .TENPERRTURE . 
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APPENDIX II 


ro.gram to Calculate Electron 
Densities and Temperatures 
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C THIS PROGRAM SOLVES A SYSTEM OF NORLIKEAR EQUATIONS EY 
C SUCCESSIVE ITERATIONS. A SYSTEM OF SEVEN SIMULTANEOUS 
C EQUATIONS HAS BEE.N REDUCED ANALYTICALLY TO TUO NONLINEAR, 

C SIMULTANEOUS EQUATIONS IN ELECTRON DENSITY AND ELECTRON 
C temperature. in EACH ITERATION STEP THE SSP SUBROUTINE 
C RTHI, SUPPLIED BY IBM IS USED TO SOLVE A NONLINEAR EQUATION*. 

C THE SUBROUTINE RTMI SOLVES FOR A ROOT OF AN EXTERNAL FUNCTION 
C SET EQUAL TO ZERO. 

COMMON DE(3) ,XKEO,TEO,XAB(3) , PH , PA , P , V , S IG 1 3 , FLX , FF , SIG 1 ,B1 1 ,B2 ,GA 
"M,S1G4 • 

• REAL YZ(610,9) 

EXTERNAL FCT 
EXTERNAL FCN 

C FORMAT PRINTED OUTPUT 

500 PORMATC' ' ,3X,'P(HE)=' ,Ei4.7,3X,"P(AR)=' , E 1 4 . 7 , 3X , ' P( TOTAL) = ' ,E14. 
*-7) - . 

550 FORMAT (' METASTABLE DENSITY = E 14,7) 

570 FORMAT (' ELECTRON DENSITY = El 4. 7) 

59 0 FORHAT(' ELECTRON TEMPERATURE= ' , El 4 . 7 ) 

610 FORLEAT(' AR2 = ',E14.7) 

611 FORMAT(' AR=',E14.7) 

612 FORMAT(' REC0MBINAT10N=' , E14 . 7 ) 

613 FORMATC' HE2=',E14.7) 

6-14 ' FORMATC' HE-',E14.7) 

620 FORMATC' ','***** NO CONVERGENCE *****') 

640 FORMATC' GARBAGE GARBAGE GARBAGE GARBAGE')^ 

670 FORMATC '1') 

680 FORMAT ('0') 

C DEFINE PARAMETERS AND INITIALIZE VARIABLES 

■N2=l 
N=2 

XNE0=1 .0E12- 
TEO=305. 

PH-1504.8 
X=-.6l ■ 

PA-X*PH 


P-PA-FPll 
PAl-PA 
PHI -PH 
Pl-P 
V-6.0E4 



■ SIG13-5.0E-21 
FLX-1 .E17 
FF-2.3E4 
SIGl-lE-17 
B1 1=1 .8E-9 
B2-2.2E-10 
GAM-2 .4E-29 
SIG4-7.0E-16 

C THE FOLLOWING SOLVES AN ENERGY' PARTITION EQUATION FOP^ THE 

C ENERGY TRANSFERRED TO THE KE AP-N AXWELLI AN PART OF A NON- 
C MAXWELLIAN ELECTRON ENERGY DISTRIBUTION FUNCTION 
200 AA-6 .93E“5*XNE0 


P-PH+PA 
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33=1 ,34E'-1 1*3 I55E1 9*P 
CC=2.36E-11*3»55E19*P 
,AAA=SQRT(AA/BB) 

SK « < AA/CC)**( 1 . /3 . ) 

A 3 C= A A A* ( A TA W ( 3 / A AA ) -AT AN ( 1 . 2 9 E- 4 * TEO?* A A A ) ) 

A B B = . 6 7 * ( X K * * 2 ) * ( ( . 5 * A L 0 G ( 1 - C 3 * X K * 1 . 7 3 / ( ( >; K+ 1 . 7 3 ) 
->C( (3.46-XK)/(1.73*XK)))} 

H=0 


"2))))+(l .73*ATA 


SN=0 

250 DO 295 K=10,20,5 
KN=0 

lF(l-(3*XK*SQRT(FLOAT(N))/(XX-!-SQRr(FLOAT(N))**2)))290,29 0»285 
. 285 XABCKN) = .67*(XK**2)*(.5*AL0G(1-(3*XK*SQRTCFL0ATCW))/ (XK+SQRTCFI-OAT 
T(N))**2)) ) + (l .73*ATAN( (<2*SQRT(FL0AT(N)) )-XK)/ (XK*1 .73)))) 

BE<H)=ABC+XAB(NN)“ABB 

IFLAG“0 

GO TO 295 

■ 290 K=M+1 

DE<H)=ABC 
IFLAG=1 - 
295 CONTINUE 

C CALCULATE ELECTRON TEMPERATURE 
S=FCT(TEO) 

IF (S) 1011,1012,1013 

C CALL SUBROUTINE TO SEARCH FOR RIGHT-HAND BOUNDARY USED IN RTMI 

1011 CALL TP0S(S ,S1 ,XKI ,XLI) 

GO TO 36 

■ 1012 TEN=TEO ' 

GO TO 37 

C CALL SUBROUTINE TO SEARCH FOP. LEFT-HAND BOUNDARY USED IK RTNl 

■; 1013 CALL TNEG(S,S1 ,XRI,XLI) 

C CALL SUBROUTINE TO SOLVE ENERGY BALANCE EQUATION 

•36 CALL RTMI(XEN,F,FCX,XLI,XRI,1. ,1000,IER) 

FK=F ■ ■ ■ • ■ 

IERN=IES. 

IF .(IER-2)37, 1999,37 

C CHECK CONVERGENCE OF ENERGY BALANCE’ EQUATION 

37. IF((ABS(TEO-TEK)/TEO)-0. 000 1)7 10,710,700 
. 700 TE0=TEN 
■ . GO TO 200 

C CALCULATE ELECTRON DENSITY 
710 S=FCN(XWE0) 

IF(S)1001 ,1002,1003 

C CALL SUBROUTINE TO SEARCH FOR RIGHT-HAND BOUNDARY USED TN HTUl 

i06l CALL NP0S(S,S1 ,XRI,XLI) 

GO' TO. 56 

1002 xnen-xkeo 
CO TO 57 

C CALL SUBROUTINE TO SEARCH FOR LEFT-HAND BOUNDARY USED' IK RTHI 

1003 CALL NKEG(S,S1 ,XRI,XLI) 

C CALL SUBROUTINE TO SOLVE ELECTRON DENSITY CONTINUITY EQUATION 

56 CALL RTN1(XNEK,F,FCK,XL1,XRI,1 .,1000,IER) 
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IF (1ERK~2)57, 1999,57 

C CHECK CONVERGENCE OF ELECTRON DENSITY CONTINUITY EQUATION 

57 IF((ABS(XNEO-XNEN)/XNEO)rO. 0001)465, 465, 460 
460 XNEO=XNEK 
GO TO 200 
1999 L'RITE(6,620) 

465 IF ( IFLAG)480,480,470 
470 WRITE(6,640) 

CALCULATE lOK AND KETASTABLE DENSITIES FROM SOLUTIONS FOR 
ELECTRON TEMPERATURE AND ELECTRON DENSITY 
480 ALF9=4.0E-20*(300,/TEN)**4*XNEN-M . 8E-10* ( 300 , /TEN)*P 
ALFD=9^ 1E"7*(300./TEN)**.61 
ALF12=7.0E-1 1*SQRT(TEN) 

C1=SIG13*FLX*3. 55EI6*PH*FF/(7 8.*PK**2+SIG1*V*3. 55E16*PA) 
C2=.5*B11/(78.*PH**2-PSIG1*V*3 .55E16*PA) 

C3=78*PH**2/(B2*3-55E16*PA+GAM*1 ,26025E33*PA*PH+ALF9*XNEN) 

C4 = (SIGl*V*3.55Ei6*PA+(B2-6GAM*3.55E16*PH)*3,55E16*P7i*C3)/(ALF9*XNE 
“M+376.*PH*PA) 

C5=SIG4*V*3.55E16*PA/(ALF9*XNEK+376.*PA*PH) 

AA10=B1 1-(27.3*ALF9^PH*PH*B11*XNEN/(7S*PK*PH+3.55E~1*V*PA)/ (ALF9*X 
-KEN+2.2E-10*3 .55E16*PA-P2.4E-29*i .26025E3 3*PH*PA) ) 

BB10=^1 .491E6*PA-PALF1 2*XNEN-5-. 6*PE*PH 

CC10=-SIG13*FLX*3.55E16*PH*FF/,56~C54.6*ALF9*PH*PH*SIG13*FLX*3.55E 
~i6*PH*FF*XNEN/(7 8*PH*PH+3>55E-'l*V*PA')/(ALF9*XHEN+2.2E-10*3.55El6*P 
'-A-P2.4E-29*! .260 25E33*PA*PH)) 
C10=(-BB10-fSQRT(BB10**2-4*AA10*CC10))/2/AA10 ' 

490 HE=Cl-!-C2*C10*C10 

HE2=C3*HE 
AR=C4*HE+C5*C10 

FIHIS=ALF9*AR*XNEN/(2.,E“9*3.55E1 6*PA+( ( ALF9*XKEN) ) ) 

AR2=(2E-9*FIKIS*3*55E16*PA+376.*PH*PA*AR)/(ALFD*XKEN) 

XKE=AR4-AR2-J-HE+HE2 

1F( ABS (XNE-XNEN) /XNEK- .01 )15,15,10 

10 PH=PKI-(HE-PHE2)/.3.55E16 

PA=PA1-CAR+AR:2) /3‘. 55E1 6 
•XKE=XKEN 
GO TO. 200 

C‘ PRINT RESULTS 
15 1IRITEC6,680). 

PRITE(6,500)PH1 ,PA1 ,P1 

URITE(6,550)C10 

WKITE(6 ,570)XKEK 

L'RITE(6,590)TEN 

URITEC6,611) AR 

,yRITE(6 ,610) AR2 

VJRITE(6,613)1IE2 

1;RITE(6 ,614)HE 

V/RITEC6,570)XKE 

STOP 

END 


C 


SUBROUTINE TPOS C S , SI ,XK I , XLl) 

THIS SUBROUTINE SEARCHES FOR THE RIGHT-HAKD BOUNDARY USED 
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c 

c 

c 

c 

c 


BY RTHl IN SOLVING THE ENERGY BALANCE EOUQTION. 

THE SEARCI: IS PERFORMED BY INCREMENTING THE ARGUMEMT OF 
THE EXTERNAL FUNCTION FCT IN THE DIRECTION OF INCREASING VALUE 
OF THE FUNCTION UNTIL A POSITIVE VALUE IS ENCOUKTEKED. C0NTR03 
IS THEN returned TO THE MAIN PROGRAL'. 


COMMON DE(3) > XNEO > TEO , XAB ( 3) 
-M,SIGA 
11=0 
J1 = 0 


PH,PA,P>V,S1G13^FLX, 


FF,SIG1,B11 »B2,GA 


• TE=1.01*TEO 


TE1=TE0 

S=FCT(TEO) 

S1=FCT(TE) 

IF (SI) 1004, 1004, 1028 . 

1004 IF (S~S1)1007, 1005, 1005 
997 TE1=TE0 

1005 TE=TE1/1^01 
S=FCT(TE) 

IF(S)1001, 1001, 1028 
1001 I1 = I1-M 

IF (Il“800)1006, 1006 ,1003 
1003 IF (31-800)999,999,1028 

1006 TE1=TE 

GO TO 1005 
999 TE1=TE0 


1007 TE=TE1*1.01-. 

S=FCT(TE) 

IF(S) 1002, 1002,1028 
,1002 J1=J1-M 

I F(J1~800) 1008,1008 , 1009 
1009 IF (11-800)997,997,1028 
1008 TE1=TE 

GO TO 1007 
1028 XRI=TE 

■ >:li=te& 


RETURK 

END- 


SUBROUTINE TKEG( S , SI ,XRI , XLI) 

C THIS SUBROUTINE SEARCHES FOR THE LEFT-HAND BOUNDARY USED 

C BY RTMI IN SOLVING THE ENERGY BALANCE EQUQTION. ’ 

C THE SEARCH IS PERFORMED BY IKCREMENTIKG THE ARGUMENT OF 

C THE EXTERNAL FUNCTION FCT IK THE DIRECTION OF DECREASING VALUE 
C OF THE FUNCTION UNTIL A NEGATIVE VALUE IS ENCOUNTERED. CONTROL 
C IS THEN RETURNED TO THE MAIN PROGRAM. 

COMMON DEC 3 ) ,XNEO, TEO , XAB ( 3 ) ,PH, PA ,P , V , SIGl 3 , FLX , FF , SIGl , B1 1 , B2 , GA 
-M,SIG4 
.31 = 0 
11 = 0 ■ 

TE1=TE0 
TE=1 .01*TEO 
S=FCT(TEO) 

S1=FCT(TE) 
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IF (Sl)1027»999,999 

999 IF (S-Sl)lOOO^IOOO, 1005 
998 TEl^^TEO 

1000 TE=TE1/ 1 . 01 ■ 

S=FCT(TE) 

IF (5)1027,1007,1007 
1007 11=11+1 

IF (il'-800)1001 , 1001 ,1002 

1002 IF (Jl“800)997, 997, 1027 

1001 TE1=TE 

GO TO 1000 
997 TE1=TE0 
100-5 TE=TE1*1.01 
S=FCT(TE) 

IF(S) 1027,1003,1003 

1003 Jl=Jl+i 

IF (Jl“800)1006, 1006,1004 

1004 IF (li'-800)99’8, 998, 1027 
1006 TE1=TE 

GO TO 1005 
1027 XLI=TEO 
XRI=TE 
.RETURN 
END 


SUBROUTINE NP0S(S,’S1 ,XRI,XLI) 

C. THIS SUBROUTINE SEARCHES FOR THE RIGHT-HAND BOUNDARY USED 

C BY RTl-II IN SOLVING THE ELECTRON DENSITY CONTINUITY EQUQTION^ 

C ■ .THE SEARCH IS PERFORMED BY INCREMENTING THE ARGUMENT 01- 
C THE EXTERNAL FUNCTION FCK IN TEE DIRECTION OF INCREASING VALUE 
C OF THE FUNCTION UNTIL A POSITIVE VALUE IS ENCOUNTERED. CONTROL 
C- IS. THEN RETURNED TO THE MAIN PROGRAM. 

COMMON DEC3) ,XKE0,TE0,XAB(3) , PK , PA , P , V , SIGl 3 , FLX, FF , SIGl , El 1 , B2 
-M,SIG4 
•JI=0 

ii=o 

XNE^1.01*XKEO 

XHE1=XNE0 

S=FCN(XNEO) 

S1=FCK(XKE) 

IF (31)11,11,28 

11 IF (S-Si)17,12,12 

7 XRE1=XNE0 

12 XNE=XNE1/1 .01 
S=FCK(XNE) ' 

IF(S)10,10,28 

10 11 = 11+1 

- IF (11-1000)13,13,8 

8 1F(J1-1000)6,6,28 

13 XNE1=XNE 
GO TO 12 

6 XNE1=XNE0 
17 XNE=XNE1*1 .01 



58 


S=FCN(XKE) 

IF(S)16> 16,28 
16 = 

IF (Jl-1000)18,18,19 
19 IF(Il-lOOO)? ,7 ,28 
IS XKE1=XNE 
CO TO 17 
28 XKI=XNE 
XLI-XNEO 
RETURN 
END 


SUBROUTINE KNEG( S , S 1 ,XRI , XLI) 

C THIS SUBROUTINE SEARCHES FOR THE LEFT-HAND BOUNDARY USED 

C BY RTHI IN SOLVING THE ELECTRON DENSITY CONTINUITY EQUQTION, 

C THE SEARCH- IS PEEFORHED BY INCREMENTING THE ARGUMENT OF 

C THE EXTERNAL FUNCTION FCT IK THE DIRECTION OF DECREASING VALUE 
C OF THE FUNCTION UNTIL A NEGATIVE VALUE IS ENCOUNTERED. CONTROL 
C IS THEN RETURNED TO THE MAIN PROGRAM. 

COMMON DE(3) ,XNEO,TEO,XAB(3) ,PH , PA ,P , V , S IGl 3 , FLX , FF , SIGl , B1 i ,B2 , G A 
-li,SIG4 
J1=0 
11=0 

XNE=1.01*XNE0 
XNE1=XKE0 
S=FCK(XKEO) . 

S1=FCN(XNE) 

IF (Si)27,9,9 

9 IF (S-Sl)IO,'lO,15 
7 XNE1=XNE0 

10 XNE=XKE1/1 .01 
S=FCK(XNE) 

IF .(S)27,l4,14 

14 ' 11=11+1 

IF Cll-20OO)li ,ii , IS- 
IS IFCJ1-1000)8^8,27 

11 XKE1=XKE 

12 GO TO 10 
a XKEr=XNEO 

15 XKE=XNEl*i .01 
S=FCK(XKE) 

IF (S)27,IS,18 
18 J1=J1+1 

IF (J1-1000)I6, 16,17 
17 • IFCI1~1000)7,7,27 

16 XKE1=XNE 
GO TO 15 

27 XR'I=XNE 
XLI=Xlv’EO 
P.ETUP.K 
END 



n n ■ , on 


59 


FUNCTION FCT(TEN) 

EXTERNAL FUNCTION USED IlY RTMl IN SOLVING EKPIRCY BALANCE 
EQUATION 

COMNON DEC 3 ) > XKEO, TEO , XAB ( 3 ), PH , PA V , S XG 1 3 , FLX , FF , SIC 1 , B1 1 , B2 , GA 
~M,SIC4 

ALFD= 9 . 1 E-7 * ( 3 00 ► / TEN) * * . 6 1 

ALF9=4.0E-20*(300,/TEN)**4*XKEO-M . 8E~ 1 0* ( 300 . / TEN)*? 

ALF12=7.0E-11*SQRTCTEK) 

AA10=B1 1-(27.3*ALF9*PH*PP>B11*XKE0/(7 8*PU*PHF3.55E-1*V*PA)/(ALF9*X 
-NEO+2.2E-10*3/55Ei6*PA+2.4E~29*l > 26025E33*PH*PA) ) 

EB10=1 ,491E6*PA+ALF12*XNE0+>6*PH*PH 

CCIO— SIG13*FLX*3.55E16*PI1*FF/.5 6-(54.6*ALF9*PH*P11*SIG13*FLX*3.55E 
-lG*PK*FF*XNEO/(78*PH*PH+3 .55E-l*V*PA)/(ALF9*XNE0+2-2E“10*3. 5 5EI6*P 

CCl0=’~SIGi3*FLX*3.55El6*PP*FF/.5 6'-C5 4.6*ALF9*PE*Pn*SIG13*FLX*3.55E 
"A-P2.4E-29*! .26025E3 3*PA*?H)) 

Ci0=('-BB10+SQRT(BB10**2-4*AA10*CC10))/2/AAl0 

XXX=SIGi3*FLX*FF*3.55E16*PH*DE(i)“1.74E"17*ABS(TEN-300.)*SQRT(T£K) 

-*3.55E16*P*XNEOF9E-10*C10**2*DE(2)+7E~11*DE(3)*C10*XKEO*SQET(TEN) 

FCT=XXX/1 *EI5 

RETUPvN 

END 


FUNCTION FCKCXNEN) 

EXTERNAL FUNCTION USED BY RTNI IK SOLVING ELECTRON DENSITY 
CONTINUITY EQUATION 

COllKON DE(3) ,XKEO,TEO,XAB(3) >PH,PA,?,V, SIG13,FLX,FF,SIG1 ,B11,B2,GA 
“M,SIG4 

ALFD=9.1E-7*(300./TEO)**.61 

ALF9-4.0E“20*(300./TEO)**4*ZNEK+1 .8E-10*(300./T£0)*P 
ALFi2=7 .OE-il*SQRT(TEO) 

Cl-SIGi3*FLX*3 .55E16*PH*FF/(78.*PH**2+SIG1*V*3.55E16*PA) 

C2 = .5*B11/(7 8.*PH**2'I-SIG1*V*3.55E]6*?A) 

C3 = 78*PK**2/(B2*3»5 5E1 6*PA+-GA1!*T . 2 6G25E3 3*PA*PH-rALF9*XNEK) 

• e4 = CSIGl*V*-3.55E16*Px\+(B2+GAN*3.55E16*PH)*3.55E16-*PA*C3)/CALF9*XNE 
-N+37 6V*PH*PA> ■ . . . 

C5=SIG4*V=i^3'.55E16*PA/(ALF9*XNEN-F37 6.*PA*PH) • . 

AA10=BI I-C27^3*ALF9*PH*PH*B11*XNEN/(78*PH*PH+3.55E~1*V*PA)/(ALF9*X 
-NEN+2.2E-10*3.55E16*PA-P2.4E-29*1 .26025E33*PH*PA)) ' ’ ' 

BB10=1.491E6 * P A-FALF 1 2* XNExN+ . 6* PH* P II 

. CC10=-SIG13*FLX*3.55E16*PH*FF/.56-(54.6*ALF9*PH*PH*SIG13*FLX*3.55E 
-16*PH*FF*XKEN/(78*PH*PH+3.55E-1*V*PA)/(ALF9'*XNEN+2.2E~10*3.55E16*P 
-A+2.4E-29*! ,26025E33*PA*PH)) 
ClO=C~BB10-fSORT(BBlO**2-4*AA10*CC10))/2/AAlO 
'-A+2.4E-29*1.26025E33*PA*PH)) 

HEI=Ci+C2*C10*C10 
HEII=C3*HEI 
ARI=C4*HEI-FC5*Ci 0 

ARX=ALF9* ARI*XNEN/ ( 2 * F.~9*3 ► 55E 1 6*PA-F( ( ALF9*XNEN) ) ) 

A RI 1= C 2 * E- 9* ARX*3 ► 5 5E 1 6 * p a+3 7 6 . * P H* P A* AR I ) / ( ALFD*XNEN ) 

YYY=liBI-HIEII-FARI-I-ARII"XNEK 

FCN=YYY/1 .E15 

RETURN 

END 


p>A@i m 
Of 'p6or'qualsiy 
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